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Abstract 

We analyze various aspects of the recently proposed holographic theories with gen- 
eral dynamical critical exponent z and hyperscaling violation exponent 9. We first 
find the basic constraints on z, 9 from the gravity side, and compute the stress-energy 
tensor expectation values and scalar two-point functions. Massive correlators exhibit 
a nontrivial exponential behavior at long distances, controlled by 9. At short distance, 
the two-point functions become power-law, with a universal form for 9 > 0. Next, 
the calculation of the holographic entanglement entropy reveals the existence of novel 
phases which violate the area law. The entropy in these phases has a behavior that 
interpolates between that of a Fermi surface and that exhibited by systems with ex- 
tensive entanglement entropy. Finally, we describe microscopic embeddings of some 
9 metrics into full string theory models - these metrics characterize large regions 
of the parameter space of Dp-brane metrics for p ^ 3. For instance, the theory of N 
D2-branes in IIA supergravity has z = 1 and 9 = —1/3 over a wide range of scales, at 
large g s N. 
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1 Introduction 

Holography [1] is a powerful tool to study strongly interacting large N quantum field theories 
[21 El IH E]- In the holographic context, a d (spatial) dimensional quantum field theory is 
mapped to a higher-dimensional (usually (d + 2)-dimensional) gravitational theory, with the 
(d + 1) "field theory dimensions" arising as the boundary of the space-time. While the initial 
interest in concrete examples centered on applications to AdS gravity theories and their 
conformal field theory duals, the class of metrics of interest in gauge/gravity duality has been 
considerably enlarged in recent years. One simple generalisation is to consider metrics which 
can be dual to scale-invariant field theories which are, however, not conformally invariant, 
but instead enjoy a dynamical critical exponent z ^ 1 (with z = 1 reducing to the case of 
the AdS metric): 

ds 2 = -±- z dt 2 + ±(dr 2 + dx 2 ). (1.1) 

These metrics are invariant under the scaling 

t -> X z t, x-i -> Xxi, r — y Xr . (1.2) 

They arise as exact solutions of simple gravity theories coupled to appropriate matter jSJ [7] , 
with the simplest such theory also including an abelian gauge field in the bulk. This simple 
generalisation of AdS is motivated by consideration of gravity toy models of condensed matter 
systems (where Lorentz invariance needn't emerge in the deep infrared, and e.g. doping with 
charge density can naturally lead to z =^ 1)F] Such metrics have also been found as solutions 
in string theory, and supergravities which arise simply from string constructions, in [9]. 

More recently, it has been realized that by studying systems including a scalar "dilaton" 
in the bulk, one can find even larger classes of scaling metrics. Such theories have been 
studied in e.g. [101 HH H21 H31 HH UHl HH H71 HE1 UHl EOl EI] (with very similar metrics also 
characterizing the "electron stars" studied in e.g. [22J). By including both an abelian gauge 
field and the scalar dilaton, one can in particular engineer the full class of metrics [13] 



ds 2 d+2 = r - 2 ^ d (-r- 2 ^dt 2 + dr 2 + dx 2 ) . (1.3) 

These exhibit both a dynamical critical exponent z and a "hyperscaling violation" exponent 
6 [24J, as emphasized in [20J. This metric is not scale invariant, but transforms as 

ds -> X 9/d ds (1.4) 



under the scale-transformation (1.2). Very roughly speaking, in a theory with hyperscaling 



violation, the thermodynamic behaviour is as if the theory enjoyed dynamical exponent z but 

1 An alternative class of non-relativistic scaling metrics characterizing the so-called Schrodinger space- 
times, which may also have condensed matter applications, was discussed in [5j. 
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lived in d—6 dimensions; dimensional analysis is restored because such theories typically have 
a dimensionful scale that does not decouple in the infrared, and below which such behaviour 
emerges. One can then use appropriate powers of this scale - which will be denoted by r> 
- to restore naive dimensional analysis. As emphasized in [20], building on the stimulating 
explorations in the case 9 = d— 1 is a promising gravitational representation of a theory 
with a Fermi surface in terms of its leading large N thermodynamic behaviour. In this 
example, the relevant dimensionful scale is of course the Fermi momentum. 

In this paper we characterize strongly coupled quantum field theories with hyperscaling 



violation using holography. In general, the metric (1.3) may not be a good description above 
the dynamical scale rp^\ For this reason, in this work we will not assume (1.3) all the way 
to the boundary, but instead we will follow an 'effective' holographic approach in which the 
dual theory lives on a finite r slice. This is similar to an effective field theory analysis in the 
dual QFT. This has been put on a firmer footing for asymptotically AdS spacetimes in [23J. 

First, we discuss the most basic holographic features of this class of metrics: the con- 
straints on (z, 6) that follow from energy conditions in the bulk, the behavior of propagators 
for bulk fields (and the consequent behavior of correlation functions of scalar operators in the 
dual field theories), and the behavior of the stress-energy tensor. Our analysis reveals intrigu- 
ing properties of correlation functions in these theories. In a semiclassical approximation, a 
massive scalar has a correlation function of the form 

G(Ax) ~ exp (-eg j\Ax\ 9/d ^) (1.5) 

at spacelike separations (where eg > is a constant). We note the nontrivial lAxj 61 ^ depen- 
dence, as compared to a weakly coupled massive propagator, G(A) ~ exp(— m\ Ax\). On the 
other hand, away from the semiclassical limit (i.e. at small masses/short distances), there is 
a cross-over to a power-law behavior, and the propagator becomes^] 

G (^) ~ lA^d+D-g • (L6) 

That is, there is a universal ^-dependent power law, independent of m. 

In another direction, we systematically study the entanglement entropy properties of 
the dual field theories. In recent years, studies of entanglement entropy have come to the 
fore as a new technique for understanding and perhaps classifying novel phases of quantum 
field theory |25j. In general, the entanglement between a region A and its complement in a 
quantum field theory in its ground state in d spatial dimensions, is expected to scale as the 
area of dA, the boundary of the region (with a precise proportionality constant dependent 
on the UV-cutoff of the field theory) [26]. For the UV-dependent contribution, this scaling 
simply follows from locality, and has come to be known as the "area law." However, several 



For instance, in systems controlled by a UV fixed point, the metric should become scale invariant for 
r <C rp. Similarly, the hyperscaling violation regime may also be changed in the deep IR if there is an 
attractive fixed point. We will find string theory examples of these situations. 
3 This assumes 9 > 0; see <|3]for more details. 



3 



states which violate the area law have also been discussed. These include d — 1 conformal 
field theories [27], conventional Fermi liquids, which can exhibit logarithmic violation of the 
area law [28J [29] , and certain proposed non- Fermi liquid ground states [30] . In these systems, 
the area law - violating terms are not cut-off dependent. More generally, sub- leading but cut- 
off independent terms in the entanglement entropy have proved to be of significant interest 
- for instance, they can distinguish between states with different topological orders [3TJ [32] . 

Ryu and Takayanagi [33] proposed that the entanglement entropy between a region A 
and its complement in the boundary field theory, can be computed in gravity by finding the 
area of the minimal surface ending on dA and extending into the bulk geometry (in Planck 
units). While this proposal is as yet unproven, it has passed many non-trivial tests, and is 
supported by an impressive amount of circumstantial evidence. Here, we systematically study 



the entanglement properties of the class of metrics (1.3), over the full range of parameters 
z, 9 where they seem to emerge as solutions of a reasonable matter + gravity theory (i.e., one 
where the required stress-energy sources satisfy reasonable energy conditions). Entanglement 
properties of subsets of these theories were studied in [T91 [201 EI]> an d also in [33], which 
emphasized the importance of the cross-over between the area law at T = and the thermal 
entropy. One of the surprises we'll encounter is the existence of a class of theories which 
violate the area law and have universal terms in the ground-state holographic entanglement 
entropy that scale parametrically faster than the area of OA (while scaling more slowly than 
the extensive entanglement entropy expected in a theory with extensive ground-state entropy 



As a third focus, we also discuss the way some 9 ^ metrics arise in a UV-complete 
theory - string theory. Existing embeddings have been in phenomenological theories of 
Einstein/Maxwell/dilaton gravity, which are clearly applicable only over some range of scales 
(as the dilaton is varying, leading one to suspect that the description breaks down both in 
the deep IR and the far UV). Here, we simply point out that over a wide range of scales, 



the dilatonic Dp-branes (those with p ^ 3) give rise to metrics of the form (1.3) with z = 1 
but 9 ^ 00 For instance, the large N D2-brane theory, in the IIA supergravity regime, 
has 9 = —1/3. The string embedding, together with our knowledge of the properties of 
Dp-branes, provides a complete understanding of what happens in the far UV and deep IR 
regions of the phase diagram where "bottom up" descriptions break down 



2 Holographic theories with hyperscaling violation 

We begin by analyzing basic properties of theories with hyperscaling violation using holo- 
graphic techniques. In this first step, our goal will be to determine two-point functions and 
the expectation value of the energy-momentum tensor for these field theories. In the fol- 
lowing sections we will construct other observables, such as the entanglement entropy, and 
study finite temperature effects. 

4 After submitting this work, we learned that a discussion which overlaps with our own discussion of 
non-conformal D-branes appears in |14j . 
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2.1 Metrics with scale covariance 



As we reviewed before, the gravity side is characterized by a metric of the form 

ds 2 d+2 = r ~ 2 ^' d (-r-^-^dt 2 + dr 2 + dx 2 ) . (2.1) 

This is the most general metric that is spatially homogeneous and covariant under the scale 
transformations 

Xi -> Xxi , t ->■ X z t , r -> Ar , ds ->■ A e/d cfe . (2.2) 

Thus, z plays the role of the dynamical exponent, and 9 is the hyperscaling violation expo- 
nent. 

The dual (d+ l)-dimensional field theory lives on a background spacetime identified with 



a surface of constant r in (2.1). The radial coordinate is related to the energy scale of the 
dual theory. For example, an object of fixed proper energy E pr and momentum p pr redshifts 
according to 

E ( r ) = Z^eJd E Pr > P( r ) = ZJZ^Ppr • (2.3) 



When 9 < dz and 9 < d, r — > (the boundary of (2.1)) describes the UV of the dual 
QFT. Clearly, different probes give different energy-radius relations, as in AdS/CFT [37] . 
For instance, a string of fixed tension in the (d + 2) -dimensional theory has E oc \/ r z - 2e / d . 
Probe scalar fields will be discussed in $3l 



Before proceeding, it is important to point out that the metric (2.1) will only give a 
good description of the dual theory in a certain range of r, and there could be important 
corrections for r — > or very large r. Outside the range with hyperscaling violation, but 
assuming spatial and time translation symmetries and spatial rotation invariance, the metric 
will be of the more general form 

ds 2 d+2 = e 2A ^ (-e 2B ^dt 2 + dr 2 + dx 2 ) . (2.4) 

An important situation corresponds to a field theory that starts from a UV fixed point and 
then develops a scaling violation exponent 9 at long-distances. This means that the gravity 
side warp factor e 2A — > R 2 /r 2 for r — > (with R the AdS radius) and that below a cross-over 



scale tf it behaves as in (2.1). This scale then appears in the metric as an overall factor 
ds 2 oc R 2 ' /r 2 p \ and is responsible for restoring the canonical dimensions in the presence of 
hyperscaling violation]^] Finally, at scales r ^> rp the theory may flow to some other fixed 



point, develop a mass gap etc., and (2.1) would again no longer be valid. String theory 



examples that exhibit these flows will be presented in £j6j For now we will simply ignore 



these corrections and focus on the form (2.1 ), keeping in mind that it may be valid only in a 
certain window of energies. We follow an 'effective' approach where the dual theory is taken 
to live at finite r of order rp. 

In order to understand the metric properties of this class of spacetimes, notice that 



(2.1 ) is conformally equivalent to a Lifshitz geometry, as can be seen after a Weyl rescaling 
5 For instance, in models with a Fermi surface, rp is set by the Fermi momentum [19] . 
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9iiv 9[iu — Q 2 9[tu, with tt = r~ Q l d . (The scale-invariant limit is 9 = 0, which reduces 
to a Lifshitz solution.) Since a Lifshitz metric has constant curvature, the scalar curvature 
associated to (2.1 ) acquires r-dependent terms controlled by the derivative of the Weyl factor 
tt. 

The Appendix contains the Ricci and Einstein tensors for the general class of metrics 



(2.4). In particular, the Ricci tensor for the metrics (2.1) is given by 

Rtt 



(d + z-9){dz-6) 
dr 2z 

-d(d + z 2 ) + (d+ z) 



Rij 



dr 2 

(d-8)(d + z 
dr 2 



(2.5) 



The scalar curvature is then R oc r 29 / d ^ which becomes constant for 9 = as expected. 



2.2 Constraints from the null energy condition 



What types of constraints should we impose on (2.1) in order to get a physically sensible 
dual field theory? Quite generally, from the gravity side we should demand at least that the 
null energy condition (NEC) be satisfied. That is, we impose 

T^N^IT > (2.6) 



on the Einstein equations, where N^N^ = 0. Since G^ u = T^ v on-shell, from (2.5) the 
constraints from the NEC become^ 

(d-9) (d(z-l)-9) >0 

(z- l)(d + z-9) > 0. (2.7) 



The constraints (2.7) have important consequences for the allowed values of (z,6) that 
admit a consistent gravity dual. First, in a Lorentz invariant theory, z = 1 and then the 
first inequality implies that 9 < or 9 > d. Both ranges will be realized in the string theory 
constructions of ^6] On the other hand, for a scale invariant theory (9 = 0), we recover the 
known result z > 1. 

Theories with 9 = d — 1 are of interest since they give holographic realizations of theories 
with several of the properties of Fermi surfaces [T9l 120] . The NEC then requires that the 
dynamical exponent satisfies 

z>2-l/d, (2.8) 



6 For the general metric (2.4), the two independent null vectors are 

AT* = , N r = e~ A cos if , N l = e~ A sin tp, 

where tp = or tt/2. 
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in order to have a consistent gravity description. More generally, in ^4] we will find that for 
systems with 

d-l<9<d, (2.9) 

the entanglement entropy exhibits new violations of the area law. These cases can be realized 
for a dynamical exponent that satisfies z > 1 — 6/d. The limit 9 = d will correspond to an 
extensive violation of the entanglement entropy, and requires z > 1 or z < 0. 

Notice that in theories with hyperscaling violation the NEC can be satisfied for z < 1, 
while this range of dynamical exponents is forbidden if 9 = [38J. In particular, {z < 0, 9 > 



d} gives a consistent solution to (2.7), as well as {0 < z < 1, 9 > d + z}. Notice that, 
just based on the NEC, the range 9 > d is allowed. Clearly more information is needed to 
determine whether the above choices lead to physically consistent theories - in particular 
we will argue below that 9 > d leads to instabilities in the gravity side. In what follows 
we continue this analysis using holographic techniques to calculate correlation functions, 
entanglement entropy and thermal effects. It would also be interesting to derive conditions 
for the existence of these theories directly in the field theory side. 

2.3 Massive propagators 

The next step is to calculate two-point functions (O(x)O(x')), where O is some operator in 
the dual theory. We will consider an operator that can be described by a scalar field in the 
bulk. The simplest correlation functions correspond to massive propagators in the bulk, since 
in the semiclassical approximation this is given in terms of the geodesic distance traveled by 
a particle of mass m. In AdS/CFT, massive bulk propagators give power-law CFT Green's 
functions because of the r-dependent warp factor. The geodesic distance is minimized by 
moving into the bulk, and this turns an exponential into a power-law; see e.g. |39j. 

Let us now calculate correlation functions in the semiclassical approximation, for the 



class of metric (2.1), in the range 9 < d. The full correlator away from the semiclassical 
limit will be studied in £j3j The particle geodesic describing the semiclassical trajectory is 
obtained by extremizing the action 



S = -m dXr-V-^V-r-^-Vi 2 + r 2 + x 2 (2.10) 



where A is the worldline coordinate and a 'dot' denotes a derivative with respect to A. The 
propagator between x = (t,xi) and x' = (t', on a fixed r = e surface is then given by 

G e (x,x')~e SM (2.11) 

with the conditions (x(0) = x, r(0) = e) and (x(l) = x', r(l) = e). Because of time and space 
translation invariance, the propagator only depends on At = t! — t and Axj = x\ — Xj. Here 
the cutoff e ~ r F , but otherwise it is left unspecified. 

Consider first the case of spacelike propagation, with At = 0. Choosing A = r gives 

S=-m f drr-( d - 0)/d y/l + {dx/dr) 2 . (2.12) 
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Since the momentum conjugate to x is conserved, the equation of motion is integrated to 



dx 
dr 



Here r t is the turning point for the geodesic, dr /dx 

|Ax| 



w ' =. (2.13) 

1 _ ( r / rt )2(d-0)/d v ' 

0. It is related to Ax by integrating 



(p3b: 



r=r t 



2d- 9 
2{d-9) 



r 



2(d-0) 



(2.14) 



Plugging (2.13) into (2.10), we obtain the geodesic distance 

d 



S= 2 -fme 9 / d --c e ,dm 



Ax\ 6 ' d 



(2.15) 



where we have neglected higher powers of e and defined 



ce,d = 




(d-9)/d 



(2.16) 



Thus, the propagator in the semiclassical approximation becomes 



G e (Ax) ~ exp 



2d 
~9 



m e 



e/d 



The approximation holds in the regime 

m \Ax\ 



exp 



> 1 



Cg d m \ Ax\ 



(2.17) 



(2.18) 



in units of the cross-over scale rp. 

As a check, in the scale-invariant limit 8 = 0, the integrals for r t and S give logarithms 
instead of powers, and we should replace e e ^ d /(9/d) — > loge (and similarly for r t ). Then 



G t (Ax) ~ exp 



m log 



Ax I 



I Ax I 



(2.19) 



reproducing the expected CFT power-law behavior. 

The correlator (2.17) reveals interesting properties about the dual field theory with hy- 
perscaling violation in the WKB regime. First, it has an exponential (rather than power-law) 
dependence on |Ax|, showing that the dual theory is massive, with a nontrivial RG evolu- 
tion - at least for operators dual to massive scalars. However, the usual weakly coupled 
decay ~ exp(— |Ax|) is now replaced by a nontrivial ^-dependent exponent. For 9 > the 
propagator decays exponentially at large distances with an exponent |Ax| 9//d . For 9 < the 
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propagator appears to approach a constant value at large distances, but it is outside the 
regime (2.18) where the semiclassical approximation is validj^] 

Anticipating our results of we point out that away from the semiclassical limit (or 
more generally for operators dual to massless scalars) the propagators do exhibit power-law 
behavior, with a power that includes a shift by 9. 

So far our calculations have been for a spacelike geodesic; similar computations lead to 
the correlator for a timelike path, now with nontrivial z dependence. Working in Euclidean 
time, r = it, the value of the action becomes, 



S 



2d 



-me 



e/d 



m 



dz 

T 



Cd,e„ 



\AT\ 0/dz 



where 



[z-9/d)T 



Cdfi,, 



dz-1 



2{z-0/d) 



e/d 



(2.20) 



(2.21) 



2(z-d/d) / 

This is valid in the range < 6/d < z. We see that for the case z = 1, this reduces to the 
solution for the spacelike geodesic. The propagator for a timelike path is 



G t (Ar) ~ exp 



2d 



m e 



exp 



dz\ 
~0 



e/dz 



Cdfi, 



m I At I 



in the regime where 



m 



Ar \e/dz > L 



(2.22) 



(2.23) 



The propagator for an arbitrary geodesic is in general a function of both \Ax\ and |At|. 
Now that we have discussed the two specific extremes, we briefly discuss the general solu- 
tion. The differential equations cannot be solved analytically for arbitrary d, 8, z, but can in 
principle be solved numerically for specific values of the critical exponents. We outline this 
procedure in appendix [B} 



2.4 Holographic energy-momentum tensor 

Another important object that characterizes the dual QFT is the expectation value of the 
energy-momentum tensor. It contains information about the number of degrees of freedom 
(e.g. the central charge in a 2d theory) and other conformal anomalies. In order to calculate 
the stress tensor, we need a method that can be applied locally to a radial slice, and which 



does not assume an asymptotic AdS structure - after all, the metric (2.1) may give a good 



description of the QFT only in an open range of radial scales. An adequate method for this 



7 In the string theory realizations of fj(3j the theories with 8 < eventually exit the regime with hyperscaling 
violation, modifying this propagator behavior. 
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case is to compute the Brown- York stress tensor [40J on a radial slice, and identify it with 
the expectation value of the energy- momentum tensor in the dual theory jHJ H2J j^] 

The basic idea is as follows. Consider a hypersurface of constant r = r c and let n M be 
the unit normal vector to this timelike surface. For us, r c ~ 0{rp). The induced metric is 

lnv = 9nv - n^nv , (2.24) 
and the extrinsic curvature is given by 

= -7/V p n„ . (2.25) 

Since r c will be taken to be finite (for instance, of order of the cross-over scale), counterterms 
and regularization issues will be ignored. Then the quasilocal stress tensor [ID] is 



T^ = K^ u -^ u K p p , (2.26) 

ignoring a dimensionful constant. 

The AdS/CFT correspondence relates the expectation value of the stress tensor (T pu ) in 
the dual theory to the limit of the quasilocal stress tensor as r c — > (the boundary): 

V^hh^(T pu ) = lim ^""V. (2-27) 

r c ->0 

where is the background QFT metric, which is related to 7^ by a conformal transfor- 



mation. Our proposal is to extend this relation to metrics of the form (2.1) at finite r, and 
use to determine {T pv ). 

A radial slice at r = r c has an induced metric 

lvbV dx»dx v = r - 2{d - e)ld {-r- 2{z - l) dt 2 + dx 2 ) , (2.28) 

and the background QFT metric is given by h^ v = r^ d 6 j/j, u - Let us first consider for 
simplicity the Lorentz invariant case. A metric of the form 

ds 2 = dw 2 + h{w)\ u dx^dx p (2.29) 

has extrinsic curvature = —h(w)d w h(w) r]^ at constant to. Applying this to our case 



and using (2.27) obtains 

(T»u) = -^r e K„. (2-30) 

c 

In the more general case of z ^ 1, we break Lorentz invariance and the nontrivial components 
of the energy-momentum tensor can be determined as 

d — 6 zd — 6 

(T 00 ) = - d+1 _ e Aioo , (T y ) = - d+1 _ e h i:j . (2.31) 

c c 



3 This has been recently applied to dc Sitter and FRW solutions, which naturally have a radial cutoff, 
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In a (d+ l)-dimensional CFT, the energy-momentum tensor is an operator of conformal 
weight d + 1, so we expect (T^) ~ h^ u /r d+1 . More precisely, a nonvanishing one-point 
function is obtained by placing the CFT on a curved background of constant curvature, and 
then (T^) ~ h^ u /R d+1 with R the curvature radius. Obtaining this from the gravity side 
requires adding counterterms and taking the limit r c — > 0. Here we are working at finite r c 
and we ignore these subtraction terms, since in general the metric with hyperscaling violation 
is not valid near the boundary. 

Hyperscaling violation has then the effect of shifting the energy-momentum tensor one- 
point function to h^ u /r d+1 ~ 9 . A similar result will be obtained in correlators of marginal 
operators below. From this point of view, a possible interpretation is that 9 reflects a nonzero 
scaling dimension for the vacuumj^] However, the effects of 6 in the field theory are probably 
more complicated than just a universal shift in the vacuum. We will return to these points 
infj 



3 Dynamics of scalar operators 

Having understood the basic properties of holographic theories with hyperscaling violation, 
in this section we study in detail operators that are described by bulk scalar fields with 
action 

S = -\ j d d+2 x y/g (<r<9 M <R0 + m 2 2 ) . (3.1) 
In particular, we will analyze two-point functions valid for arbitrary (not necessarily large) 



mass m, where the WKB approximation of £2.3 is not applicable 



3.1 Scalar field solution 



The equation of motion for a scalar field with mass m in the background (2.1) is 



Let us first consider the behavior of <fi at small r. Starting from an ansatz <fi ~ r a , we find 
that the <9j, d t , and m 2 terms are all subdominant at small r if z > and 9 > 0. In this 
case, we can solve the equation at leading order in r, which gives a = 0or<i — 6 + z. This 
means that when we impose the incoming boundary condition at r = oo (or in the Euclidean 
picture, the regularity condition), the full solution has the following expansion around r = 0: 

= 1 + .-. + G(k, u)r d - 9+z {l + •••), (3.3) 

where we have Fourier transformed in the t and x directions, and ■ ■ ■ refers to higher order 
terms in r. 



3 We thank S. Shenker for interesting remarks along these lines. 
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The behavior in (3.3) should be contrasted to the case of 9 = 0, in which the mass term 



becomes one of the leading contributions, and we are back to the Lifshitz or AdS case of 



The momentum-space two-point function on the boundary for the operator dual to is 
given by the coefficient function G(k,u) in the expansion above [IT]. We will analyze its 
behavior in the next few subsections, while solving it exactly in a few special cases. 

3.2 Massless case 

For simplicity we will consider the z = 1 case where we recover Lorentz invariance. The 



equation of motion (3.2) becomes exactly solvable for a massless scalar: 



d 2 - —d r -k 2 )<i> = 0, (3.4) 
r / 

with k = (oj, k). The solution that satisfies the correct boundary condition at r = oo is 

<P=(krf d - e+1 V 2 K {d _e + i)/2(kr) . (3.5) 



Note that we have normalized the solution at the boundary according to (3.3), up to a 
numerical factor that does not depend on k. Expanding the modified Bessel function, we 
find (again up to a ^-independent constant) 

G(k) ~ k d - e+1 . (3.6) 

Fourier transforming back to position space, we find the two-point function to be 

/rfd+l u 1 

Here O is a marginal operator dual to the massless <fi in the bulk. We find that the dimension 
of this marginal operator is shifted by 9. 

3.3 Massive case: a scaling argument 

In the more general case where the mass is nonzero, we cannot solve the scalar equation of 
motion 

9 r - — dr -k 2 - = (3.8) 

in closed form (except for special values of 9 which we will discuss in the next subsection). 
However, we note a scaling symmetry in the equation 

r^Ar, k^k/X, m^m/X e/d , (3.9) 
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under which the coefficient function G(k) should transform as 

G(k; m) = X d - e+1 G(k/X; m/\ e/d ) (3.10) 



in order to keep (3.3) invariant. 

We immediately observe that in the massless case G(k) ~ k +1 , by setting A = k. This 
agrees with our results in the previous subsection. 

We also find that for positive 9, the mass term become unimportant at short distances, 



and the UV behavior of the massive two-point function is given by the massless results (3.6 



3.7). The long-distance behavior of the massive two-point function is given by the WKB 



approximation of £2^3 We will verify these statements in a few exactly solvable cases in the 
next subsection. 

When 9 is negative, the mass term becomes unimportant at long distances, and the IR 
behavior of the massive two-point function is given by the massless results. 

We have restricted ourselves to the z = 1 case here, but our results apply more generally 
for z 7^ 1 as well. In that case the scaling symmetry is 

r-»Ar, k^k/X, u-+cu/X z , m^m/X e/d , (3.11) 

and the momentum-space two-point function transforms as 

G(k, w; m) = X d - e+z G{k/X, u/X z ; m/X e/d ) . (3.12) 

Fourier transforming back to position space, we have 

G(Ax, At; m) = X 2{d+z) - e G(XAx, X z At; m/X e/d ) . (3.13) 

The equal-time two-point function in the massless case is therefore given by 

- \^\i +z) -e ■ ( 3 - 14 ) 

3.4 Some special cases 



The equation of motion (3.8) simplifies and becomes solvable in some special cases, which 
we now discuss. 

3.4-1 The case 9 = d 

The solution that satisfies the correct normalization and boundary condition is 

<p = exp ^— y/k 2 + m 2 r^ . (3.15) 

The two-point function in momentum space is therefore 

G(k) = Vk 2 + m 2 . (3.16) 
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At short distance, the two-point function is dominated by the large-fc behavior G(k) ~ k, 
and we have 

1 



(O(x)O(x')) 



I r^tf I d~\~ t 2i 



for small \x — x'\ 



(3.17) 



which agrees with (3.7) for 9 = d. 



At long distance, the two-point function in position space 

d d+l k 



(O(x)O(x')) 



-Vk 2 + m 2 e ik < x - x,) 



(3.18) 



(27r) d+1 

can be shown to decay as ^- m \ x - x '\ by the saddle point approximation. This agrees with 
pTT7|) for 9 = d. 



3.4.2 The case 9 = d/2 

The equation of motion ( |3.8[ ) is exactly solvable for 9 = d/2 in terms of the confluent 
hypergeometric function. A special case of this kind is d — 2 and 9 = 1, which is a candidate 
holographic realization of a Fermi surface in 2 + 1 dimensions [20] . 

For 9 = d/2 with general d, the solution with the correct normalization and boundary 
condition is 

, _ r(i + i + if) e - w+iu (* + >4,- d , 2kr ) ( 3.i9) 



r(i + f ) ' V 4 2k ' 2 



2 , 



(3,0) 

from which we read off the two-point function in momentum space: 



G(k) = (2k) 1+ ^ 2 1 Z ? k ' ■ (3-21) 

r(i + 0r(-f + g) 



Again, the short- distance behavior of the two-point function is given by G(/c) ~ z^+d/ 2 a t 
large fc, and agrees with (3.7) for 9 = d/2. 



When d is even, the first gamma function in the numerator of (3.21 ) diverges and indicates 
the appearance of logarithmic terms in k. The two-point function in momentum space 
becomes 

G(k) ~ k 2 \ogk (3.22) 
at large k, which gives us G(/\x) ~ l/|Ax| 3 . 
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4 Entanglement entropy 



In this section we evaluate the entanglement entropy for systems with hyperscaling violation, 
according to the holographic proposal of Ryu and Takayanagi Our main result is the 



entropy formula (4.24) for theories with arbitrary (z,9). We will use this to probe various 
properties of these theories, including ground state degeneracies and the appearance of Fermi 
surfaces. Our study will reveal novel phases for d — 1 < 9 < d, which feature violations of 
the area law that interpolate between logarithmic and linear behaviors. 

A natural question in a holographic study of entanglement entropy is which system 
extremizes the entanglement entropy over a given class of metrics. This question is motivated 
in the following sense: one measure of strong correlation is ground-state entanglement, and 
it is well known that some of the most interesting systems (Fermi liquids, non-Fermi liquids) 
have entanglement which scales more quickly with system size than 'typical' systems. It 
is therefore worthwhile to ask, does holography indicate new phases (dual to new bulk 
metrics) with equally large or larger anomalous ground-state entanglement? This was one 
of our original motivations in this analysis. In §4.4| we address this question for metrics with 
hyperscaling violation, finding that 9 — d — 1 is the only local extremum. This implies that 
systems with a Fermi surface minimize the entanglement. 

4.1 General analysis 

Before computing the entanglement entropy in systems with hyperscaling violation, let us 
discuss the more general class of metrics 

ds 2 d+2 = e 2A ^ (-e 2B ^dt 2 + dr 2 + dx 2 ) . (4.1) 

We will first calculate the entanglement entropy across a strip. This is the simplest case to 
analyze. Then it will be argued that the same behavior is found for general entanglement 
regions when their diameter is large. 

Let us then begin by computing the entanglement entropy for a strip 

-/<£!</, 0<x t <L, i = 2, . . . , d (4.2) 

in the limit I <C L. We focus on the case of 9 < d so the strip is located on a UV slice at 
r = e. The profile of the surface in the bulk is r = r(a; 1 ), and its area is 

A=L™ r^W./i+l^i) . (1.3, 




We have inverted x\ = x\(r) to make the conserved momentum manifest, and the turning 
point r t corresponds to drjdx\\ rt =0. To obtain the entanglement entropy we need to 
extremize A and evaluate it on the dominant trajectory. 
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The calculation follows the same steps as those of £2^3 for the particle geodesic. Extrem- 
izing A and using the conserved momentum obtains 



-d{A{r)-A(r t )) 



dr y/l _ e -2d(A(r)-A(r t )) ' 

The turning point is thus fixed in terms of the length I by the integral of this expression 



(4.4) 



dr 



-d{A{r)-A(r t )) 



-2d(A(r)-A(r t )) 



(4.5) 



Finally, replacing (4.4) into (4.3) obtains the formula for the area 

p dA{r) 



A = L 



d-l 



dr 



vT 



-2d(A(r)-A(rt)) 



(4.6) 



The entanglement entropy for a strip in the general metric (2.4) is thus 



M d 



(4.7) 



with Mpi the (d + 2) -dimensional Planck constant. 
4-1.1 General entanglement regions 

While most of our analysis will be carried out explicitly for the simplest case of a strip, our 
conclusions will also apply to general entanglement surfaces. We will now establish this, by 



showing that the entanglement entropy for a general surface is given approximately by (4.6) 
both near the boundary and at long distances. 

Consider a general surface, parametrized by 



x d = <j{xi) , i = l,...,d-l 
at r = 0. The surface that extremizes the area will then be of the form 

x d = E(xi, r) , S(xi, 0) = a(xi) . 
The pullback of the bulk metric onto E gives 

dsl = e 2A ( [1 + (<9 r £) 2 ] dr 2 + 2<9 r £<9 i £ drdx { + (<% + ^S^E) dx l dx j ) 
and the area reads 



(4.8) 

(4.9) 

(4.10) 
(4.11) 
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The equation of motion implies the existence of a conserved current Jm with components 



Ji 



D dA(r) 



dA(r) 



Integrating 8m J 



M 



V / l + (d i £)2 + /(r)(<9 r £)2 

y/1 + + (9 r S) 2 ' 

over Xi, we read off the conserved charge 

as 



Qr 



d d-l xe dA(r) 



y/i + (as) 2 + (as) 



(4.12) 



(4.13) 



which generalizes the result for a strip (4.4) to an arbitrary shape. 

We will now show that (4.11 ) reduces to the case of a strip (4.6) both near the boundary 
r = e and near a 'turning point' as — > oo. First, for r — > 0, S(xj,r) may be expanded as 10 



E(xi,r) 



a(xi) + r aAxi) + ... . 



(4.14) 



The equation of motion for <j\ then requires A = 1. Then the UV part of the area is of the 
form 



A uv « (J d^x ^l + (di<r)A jdre dA ^ T) 



(4.15) 



which indeed agrees with the result for the strip ( 4.6[ ). 

Now we want to show that for regions of large area (or diameter), the long-distance part 
of the entanglement entropy also coincides with (4.6). Intuitively, when the size of the system 
is large most of the surface is deep inside the bulk, and the scaling of the entropy can be 
approximated by the behavior in the vinicity of a turning point r = r t with as — > oo. 

In more detail, we require that locally around r t , <9«S is smooth and that the combination 

e dA{r \r -r t ) 1/2 ->0 as r -> r t . 



This guarantees that Jj — > as r — > r t (see (4.12) for a definition of Jm) and hence the 
current conservation equation implies that J r ~ const near r = r t . In this approximation, 



as 



(as) 



-d(A(r)-A(r t )) 



-2d(A(r)-A( n )) 



(4.16) 



This agrees with the behavior (4.4), so the entropy from the IR region also agrees with that 
for a strip, 

e dA(r) 



A m « Q d d ~ l x y/l + (aS) 2 ^) p dr 



y/l - e -2d(A(r)-A(r t )) 



(4.17) 



Given these results, in what follows our calculations will be done explicitly for a strip, 
keeping in mind that our conclusions are valid for more general entanglement regions that 
satisfy the above criterion. 



10 A similar analysis near the boundary for a specific metric appeared recently in 
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4-1-2 Using trial surfaces 



Let us also briefly mention that the basic properties of the entanglement entropy can be 
understood by considering a simple trial surface in the bulk (e.g. a cylinder) and requiring 
that it is a stationary point. This method also applies to a general entanglement region. 

Consider a general entanglement region S defined on a slice r = e; denote its volume 
| E | and surface area by | <9£ | . We now approximate the bulk surface used in the holographic 
calculation of the entanglement entropy by a cylinder with boundary <9£, that extends from 
r = e to r = r t . The value r t is chosen such that it extremizes the entanglement entropyp] 



Starting from the general metric (4.1), the total area of this bulk cylinder is then 

A=\dZ\J t dre dA{r) + \E\e dA{rt) . (4.18) 
Requiring that r t is a stationary point gives 



A'ir.) = - (4.19) 



1 


dE\ 


d 


|E| 



where we introduced the perimeter / = |E|/|<9£|. 



Given a concrete warp factor, (4.19) determines the value of r t , and then (4.18) gives 
the approximation of the entanglement entropy by trial cylinders. The stationary point is a 
minimum or a maximum depending on the sign of A"(r): 

5 -^: = \Z\de dA ^A"{r t ), (4.20) 

evaluated at the critical point. In examples below we will find that this procedure gives a 
good qualitative understanding of the entanglement entropy. 



4.2 Entanglement entropy with hyperscaling violation 



Now we are ready to evaluate the entanglement entropy for metrics with hyperscaling vi- 
olation. It is useful to start by recalling the scale invariant 6 = case. The bulk metric 
corresponds to e 2A = R 2 /r 2 in the ansatz (4.1), with R the AdS radius. Plugging this into 



(4.5) and (4.6) obtains 



(RM, 



pi 



4(d-l) 



d-l 



r(£) 



i-l 



(4.21) 



up to higher powers of e. For d = 3, this reproduces the entanglement entropy for M 
SYM, after relating the 5d Einstein frame quantities to their lOd counterparts. 



n Of course, this approach is well-known in applications of the method of variations. See e.g. 
earlier work in the context of entanglement entropy. 



for an 
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The hyperscaling violation exponent modifies the warp factor to e 



dA(r) 



ignoring 



for now the cross-over scale). Now (4.6) gives 
A = L 



-<2-l„l-<M 
' t 



du 



U 



-(d-9) 



L 



d-1 



r 



l+d-8 

2{d-e) I r 



,l-d-\ 



,i-d+e 



2{d-6) 



1-d+e l-d+9 



and, from (4.5), the turning point is related to the length of the strip by 



■j ^ u (d~e) 
r f I du , 

lo Vl - u 2 ^ 



i+d-e 

2{d-8) 



1 



2(d-0) 



(4.22) 



(4.23) 



Restoring R and the cross-over scale rp, we thus find that the entanglement entropy 
across a strip is 

d-6 



S 



(M Pl R) d 
4(d-9- 1] 



e 



d-l 




I 



d-l 



(4.24) 



again neglecting higher powers of e. Comparing with the scale invariant answer (4.21), we 



see that the effect of the hyperscaling violation exponent is to modify the entropy by an 
additional power of (length) 61 . This can be understood directly in terms of scaling weights: 
since the metric has dimension 6/d, the entanglement entropy across a <i-dimensional region 
acquires a scaling weight 9. 

It is also useful to obtain the entanglement entropy using the method of trial surfaces 



described above. Choosing a cylinder in the bulk, (|4.18|) and (|4.19|) evaluate to 

(M Pl R) d 



trial 



4{d-6-l) 



e 
r F 



c d-l 



(d-9) 



d-d 



I 



id-i 



(4.25) 



where the trial value of the turning point is r t = (d — 9)1 and recall that / = |E|/|9E| here. 



Eq. (4.25) correctly reproduces all the physical features of (4.24); moreover when S is a strip 



both expressions exactly agree for 9 — d — 1. Eq. (4.25) of course applies to general surfaces 



indicating that our conclusions are valid beyond the case of strip-like region. 



4.3 Novel phases with d — 1 < 9 < d 



The entanglement entropy result (4.24) reveals interesting properties of the dual field theory. 



First, when 9 = d — 1 the integral in (4.22) gives a logarithmic (instead of power-law) 
dependence, so that the entropy becomes 



4 



r F I 



(4.26) 
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The cutoff can be chosen to be of order the cross-over scale, e ~ rp, where we expect the 



metric (2.1) to give a good description of the dynamics. Eq. (4.26) shows a logarithmic 



violation of the area law, signaling the appearance of a Fermi surface in the dual theory. 
This case was studied in detail by [TjJl EH] , who identified various properties of this (strongly 
coupled) Fermi surface. In particular, rp ~ kp 1 , the inverse scale of the Fermi surface. 
Another special value corresponds to 9 = d, where the metric becomes 



ds 1 



R 2 



(_ r -2(*-l) dt 2 + dr 2 + ^ 



(4.27) 



and the geometry develops an IR d factor. In this limit, the surface that bounds the entan- 
glement region does not move into the bulk, and the entropy is simply proportional to the 
volume of the entanglement region, 



S 



{MpiRf L d ' l l 



(4.28) 



This is an extensive contribution to the entanglement entropy and suggests that the dual 



theory has an extensive ground state entropy. Note that the metric (|4.27|) is not that of 
AdS2 x IR d , which shares this feature, 
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Having understood these two limits, it is clear that in the range of parameters 

d - 1 < 9 < d, 



(4.29) 



(4.24) predicts new violations of the area law that interpolate between the logarithmic and 



linear behaviors. These novel phases present various intriguing properties. To begin with, 
the entanglement entropy is finite: (4.24) does not diverge if we take the cutoff e — > 0. 



Also, in §2.2| we learned that in general these systems have a nontrivial dynamical exponent 
z > 1 — 9/d. The correlation functions computed in §^2] and [3] may also provide information 
to further characterize these phases. For instance, for massless scalars the two-point function 
is G(|Ax| ~ l/|Ax| d+3 - Q , where a = 9 - (d - 1). 

In the following sections we will derive further properties of these systems by placing 
them at finite temperature and will comment on the possible ground states that can lead to 
these new phases. 



4.4 Extremizing the entanglement entropy 



Finally, based on the result (4.24), let us determine the value of 9 that extremizes the entropy. 



We focus on the finite contribution to S in the limit where the diameter is much larger than 
rp, which is necessary in order to obtain a universal answer independent of the entanglement 
surface o 

12 Similar volume laws were discussed in [35] (and [H] for flat space holography) and related to nonlocality; 
we see no reason the metrics we are studying here are dual to non-local theories, however. 

> 0. In this case it is consistent 
also yields a universal answer. 



13 The same result is obtained if the metric ( |2.1[ ) is valid all the way to r - 
to take e ~ rp — > 0, and extremizing the first term of (4.24) with respect to 
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Extremizing the second term of (4.24) with respect to 9 shows that there is a local 
minimum at 

1 



9 C w d - 1 - - ., . . . (4.30) 
\\og(l/r F )) 

Therefore, 9 = d — 1 minimizes 5 in the limit of large diameter. This establishes that Fermi 
surfaces are local minima of the entanglement entropy]^] 

We can also ask which value of 9 gives the strongest /-dependence. Assuming 9 < d, the 
strongest rate is a linear dependence in I, when 9 = d (see also p2]). The entropy then scales 
like the volume of the entanglement region. Intuitively, this is associated to the logarithm 
of the number of degrees of freedom, which scales like the volume. 



However, just from the NEC analysis in £2.2, we found that 9 > d is allowed. In this 



range, the entanglement entropy scales faster than the volume, which is not expected to 
correspond to a QFT behavior. It is interesting to ask whether the entanglement entropy 
can reveal additional properties of this regime. For this, consider again the calculation in 



terms of trial cylinders in the bulk. At the stationary point, the entanglement entropy (4.25) 



is a good approximation to the exact answer. The "mass" for fluctuations around the critical 



point r t = (d — 9)1 can be calculated from (4.20), which now implies 



tt«— ■ 4.31 

For 9 < d, the stationary surface is a minimum; at 9 = d the bulk surface collapses into 
the r = e slice, explaining the extensive scaling. And for 9 > d the stationary point becomes 
a maximum. This suggests that gravitational backgrounds with 9 > d cannot appear as 
stable theories (at least in the regime of validity of our current analysis). ^ exhibits similar 
thermodynamic instabilities, and in the string theory realizations one finds that there is no 
well-defined decoupling limit. All these results suggest that theories with 9 > d may not be 
consistent. 



5 Thermodynamics 

So far we studied properties of QFTs with hyperscaling violation at zero temperature, such as 
correlators, the energy-momentum tensor, and entanglement entropy. We will now analyze 
finite temperature effects in these holographic theories. After obtaining the basic thermody- 
namic quantities for this class of theories, we will study the entanglement entropy at finite 
temperature and the cross-over to the thermal result. In particular, this will reveal how 
the degrees of freedom responsible for entanglement are related to those that are excited by 
thermal effects. The reader is also referred to p21 El [201 US] for related discussions. 



14 It would be interesting to study explicitly other entanglement regions. 

15 For 8 > d the UV boundary is at r — oo so one has to put the entanglement region E there. Eqs. (4.20) 



and (4.311 are still valid in this case whereas (4.19) changes sign. 
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5.1 Gravitational background at finite temperature 



Finite temperature effects are encoded in the general metric (2.4) by introducing an em- 
blackening factor fir) with 

ds 2 d+2 = e 2A ^ [~e 2B ^f{r)dt 2 + ^ + dx^ . (5.1) 

The basic property of f(r) is that it vanishes at some r = r^; the temperature is then 
proportional to a power of r^, as we explain in more detail below. In order to study finite 
temperature effects on a regime with hyperscaling violation, in the gravity side we need to 
take r F < r h . 

Starting now from the metric (2.1) with hyperscaling violation, the black hole solution 
becomes [HI [31] 

= PT (£) (~r-^f{rW + + «fa?) , (5.2) 

with 

(\ d+z-e 
k) • (5 - 3) 

Starting from a solution with f(r) = 1 and matter content general enough to allow for 



arbitrary (z,9) one can show, using the results in the Appendix, that (5.3) still gives a 



solution. Concrete examples will be presented in §61 As usual, the relation between the 
temperature and follows by expanding — r = u 2 , and demanding that near the horizon 
the metric is ds 2 ~ du 2 + u 2 dr 2 , where r = (2irT)it. The result is 



i \d + z-e\ 

Air r z h 



(5.4) 



These expressions imply that the thermal entropy, which is proportional to the area of 
the black hole, becomes 

rp(d-8)/z 

S T ~ (M Pl R) d V g — . (5.5) 



Thus, a positive specific heat imposes the condition 

d-0 



>0. (5.6) 

We see that the branch {0 < z < 1,6 > d + z} that was consistent with the NEC is 
thermodynamically unstable. On the other hand, {z < 0,9 > d} is still allowed by (5.6). It 
would be interesting to study this case in more detail to decide whether it is consistent - the 
entanglement entropy analysis of ^4] suggested an instability for all 9 > d. 



16 



This can be accomplished for instance in Einstein gravity with a scalar and gauge field. 
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Eq. (5.5) suggests that d — 9 plays the role of an effective space dimensionality for the 
dual theory. From this point of view, 9 — d — 1 yields a system living in one effective space 
dimension, i.e. a (1 + l)-dimensional theory. Recall also that for this value of 9 there is a 
logarithmic violation of the area law for the entanglement entropy. These points support 
the interpretation of 9 = d — 1 as systems with a Fermi surface [T9l [2D] . The case 9 = d 
would then correspond to a system in (0 + l)-dimensions. In ^4] we found novel phases with 
d — 1 < 9 < d that violate the area law. According to this interpretation, these would be 
systems of defects living in a fractional space dimension (between and 1). 

Notice that the behavior of the thermal entropy for 9 = d can also be obtained in systems 
with 9 < d by taking z — > oo. This is the familiar AdS2 x M. d limit of a Lifshitz metric. We 
see from the metrics that these systems are not equivalent, and they are distinguished in 
terms of field theory observables by their correlation functions. In particular, the two-point 
function for a marginal operator, 

G(x,x')= ^_ x ,^ {d+1) _ e , (5.7) 

implies that positive 9 increases the correlation between spatially separated points. There- 
fore, despite the extensive ground-state entropy, 9 = d does not in any sense correspond to 
spatially uncorrelated quantum mechanical degrees of freedom. 

5.2 Entanglement entropy and cross-over to thermal entropy 

We now study the entanglement entropy at finite temperature. This quantity is of physical 
interest since it illustrates how the degrees of freedom responsible for the entanglement 
entropy contribute to the thermal excitation. As [34J argued recently, we expect a universal 
crossover function that interpolates between the entanglement and thermal entropies. 
Finite temperature effects modify the entanglement entropy for a strip as follows: 

M d f n p dA(r) 

4 J e JK J y/l - e -2d(A(r)-A(r t )) V ; 

where r t is given in terms of the length and temperature by 

rrt -d(A(r)-A(r t )) 

1= / drf(r)- 1/2 -= =. (5.9) 

For the purpose of comparing with the thermal entropy we focus on the universal finite 



contribution to (5.8). 



Evaluating these expressions for the metrics (2.1) with hyperscaling violation obtains 



W=^Vv;+-/_(li) , Z = r,/ + (£] (5.10) 
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where we have defined the integrals 



J±(a) 



du 



u 



±(d-9) 



y/l - (au) z + d - e Vl - u 2 ( d ~ 



(5.11; 



Also, recall from (5.4) that ~ T x l 



We first need to express the turning point r t in terms of I and T and then plug this into 
«S finite- While I±(a) don't have a simple analytic form for general 9, we can analyze the 
limits of small and large temperatures explicitly. In the small temperature regime we have 
T l l z l <C 1 or, equivalently, 3> r t . The lowest order thermal correction to the entanglement 



entropy (4.24) is 



S finite OC (M Pl R) C 



d-1 



I 



1 



d-9-1 



+ ci{d-9,z)(lT 



l/z\d-6+z 



+ 



(5.12) 



with ci a positive constant that depends on d — 9 and z. So the finite contribution to 
the entanglement entropy increases by thermal effects, with a nontrivial power r p^ e+z ^ z . 
Interestingly, this dependence is in general non-analytic - for instance, the leading thermal 
correction is ~ T^ z+1 ^ z in systems with a Fermi surface (9 — d — 1). It would be interesting 
to understand the physical implications of these corrections. 

On the other hand, in the large temperature regime, — > r t , and the integrals I± are 
dominated by the pole near u — 1. Then l/rh, and expressing in terms of T 

finds 

Sfi^eOcL^ZTM/*, (5.13) 



which agrees with the thermal entropy (5.5 ). This verifies the existence of a crossover function 



that interpolates between the entanglement and thermal entropy, in theories with 9 7^ 0. 



6 String theory realizations 

In this last section we will show how some 9 7^ metrics arise from string theory. Theo- 
ries with nontrivial hyperscaling violation have been realized so far in effective actions for 
Einstein, Maxwell and dilaton fields. However, this description usually breaks down in the 
far UV or IR, so it is important to have UV completions that explain what happens at very 
short or very long distances. We will accomplish this by noting that, over a wide range of 



scales, Dp-branes with p ^ 3 give rise to metrics of the form (2.1) with z = 1 but 9 7^ 0. 
This discussion has overlap with similar remarks in [14"] . 

6.1 Black-brane solutions 

Let us first review the necessary results on black branes in ten-dimensional supergravity. 
For more details and references to the literature see e.g. [5J US]. In lOd string frame, the 
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Dp-black brane solution is 



1 5- 



ds 2 str = -4=dt* + v//^) $>, 2 + f -^ff^-d? + P 2 f-iP?-^dnl p , (6.1) 

Vf-(P) i= i j +[ *p> 

where 

The supergravity solution includes a dilaton 

e -2*W = g -2 f _( p y(p-3)/2 (63) 

and RR p-form with field strength 

/ *F P+2 = N, (6.4) 

where N is the number of D-branes. The ADM mass is M oc (8 — p)p 7 + p — p 7 S p , and 
N oc (p+p_)( 7 - p )/ 2 . 

It is convenient to introduce a new radial coordinate 

U 7- P = pi- P _ p i-P ( 6 5 ) 

and define 

p^-f = u ); v cosh 2 /3 , p 7 _- p = u 7 h - p sinh 2 /3 . (6.6) 

Here 

sinh 2 = - X - + ^ + ( Cp g s N(l s /u h y-p) 2 . (6.7) 
Then the metric and dilaton acquire the more familiar form 

ds 2 str = H{u)~W (-f(u)dt* + £ dxfj + H(u)^ (J^ + « 2 dS^ 

e^ u) = g s H(uf^l\ (6.8) 



with 

H(u) = 1 + sinh 2 ^ , f(u) = 1 - ^ ■ (6.9) 



7-p 7-p 



In order to compute the Bekenstein-Hawking entropy, we need to change to lOd Einstein 
frame. This is accomplished by rescaling 

ds 2 E = {g; l e*)-" 2 ds 2 str . (6.10) 
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The area of the horizon at u = Uh then gives an entropy 

<S bh ~ cosh f3 t\~ p , (6-11) 

in lOd Planck units. This defines a thermal entropy for the dual theory on the Dp-brane, 
where the temperature is determined by cosh/3 and u^- 

T i— . (6.12) 

cosh p Uh 

These results simplify in the limit of sm all temperature, in which case the black branes are 
nearly extremal. When Uh — >■ 0, (6.7) gives cosh 2 /? ~ sinh 2 /3 ~ l/u h ~ p . Note also that 



sinh 2 f3 u h p ~ g s N in string units. Then T ~ u h 5 p ^ 2 and 

9-p 

Sbh ~ T 5 -p . (6.13) 

The (extremal) supergravity description is valid when the curvature and dilaton are small. 
In terms of the effective 't Hooft coupling on the branes, 

u 6 p 

the solution is valid when [35J 

1 < g 2 cS < N^p . (6.15) 

At large iV and for p < 7, this gives a large range of u where the supergravity description 
can be trusted. 

Notice that the dilaton grows large for p < 2 in the deep IR, and goes outside the range 



(6.15). For example at p = 2, the theory flows into the M-theory regime. We emphasize that 



the entropy scaling (6.13) is valid when the horizon is located within the regime of validity 



of 10D supergravity/string theory (6.15), and at the corresponding range of temperatures, 



the theory exhibits hyperscaling violation. 

6.2 (p + 2)-dimensional effective theory and hyperscaling violation 

We will now compactify this theory on S 8 ~ p and show that it leads to hyperscaling violation. 
Dimensionally reducing on the sphere and changing to Einstein frame in p + 2 dimensions 
obtains 

ds 2 p+2 = u^-^ p H{u)^ (-f(u)dt 2 + j^dxj + H{u)^j . (6.16) 

This is of the general form (2.4) after the redefinition 

dr = H(u) 1/2 du. (6.17) 
Also p = d in the notation of the previous sections. 
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Taking the near horizon limit 17 of (6.16), we arrive at metrics (2.1) and (5.2) with hy- 
perscaling violation exponent 

e = P - 9 -^, (6.18) 

5 — p 

where r oc u^" 5 '' 2 . The emblackening factor f(u) also reproduces the black hole solution 



(5.2 ). As a further check, we can compute thermal effects in this effective theory and compare 



with the lOd answer. For instance, plugging (6.18) into the formula for the thermal entropy 
(5.5) indeed agrees with (6.13) for z = 1. 



We conclude that black-branes with p ^ 3 give rise to metrics with hyperscaling violation. 



This description is valid in the range of radial variables (6.15), and provides an explicit dual 
field theory realization of systems with hyperscaling violation. The field theory is given 
by SU(N) super Yang-Mills in (p + 1) dimensions, with sixteen supersymmetries; for a 
large range of energy scales and strong 't Hooft coupling, it realizes a hyperscaling violation 



exponent (6.18). Notice that 6 < for p < 4, and 9 > p for p > 6. Of course, these values 



satisfy the NEC constraints found in §2.2[ It is important to remark that the p > 6 cases do 
not, however, "decouple from gravity" and give rise to well-defined non-gravitational theories 
the way the p < 5 cases do. 

A relevant case for condensed matter applications is that of iV D2-branes, which lead to 
6 = —1/3. The string theory realization also allows us to understand the deep UV and IR 
limits, where the hyperscaling violating regime is not valid. The UV theory is given by the 
maximally supersymmetric YM theory in 2 + 1 dimensions, which is asymptotically free. In 
the IR the theory flows to a strongly coupled conformal field theory dual to AdS^ x S 7 - the 
gravity background for M2 branes. 

Besides providing an explicit string theory realization of systems with hyperscaling vio- 
lation, our results from the previous sections show interesting properties of maximally sym- 
metric Yang-Mills theories in the intermediate strongly coupled regime (6.15). For instance, 
the propagator for an operator dual to a massless scalar is 



G(x,x') 



1 



\x 



X 



/IP+2- 



9-p 
5-p 



(6.19) 



The entanglement entropy is also a probe of strong dynamics. The result (4.24) yields 

L d-i 



S 



a 



L 



d-l 



P 4 



(6.20) 



with C p a numerical constant. The entanglement entropy for D-branes was also calculated 



in [36]. Their interpretation of (6.20) was in terms of an area law and a number of degrees 



of freedom with nontrivial dependence on the RG scale, 



N pf f oc e 



(6.21; 



17 



Note that the near horizon limit u 7 p <C g s N overlaps with the supergravity range (6.15) for p < 7 
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7 Future directions 



In this paper, we have systematically analyzed the most basic holographic characteristics of 



the family of metrics (1.3), building on much earlier work. There are many directions in 
which one could imagine further developments. 

On the one hand, the classes of such metrics which are known to arise in string theory 
are still quite limited. We saw here that Dp-brane metrics in the supergravity regime, for 
p ^ 3,5, provide one class of examples. AdS2 x M 2 and Lifshitz spacetimes provide another. 
But the cases of most physical interest in conventional systems, such as the 9 = d — 1 
"Fermi-surface" - like case [20J, remain to be realized. 

On a related note, it would be interesting to interpret the 9 ^ metrics more explicitly 
in terms of dual field theories. In some of the cases with 9 < 0, it may be useful to think of 
these metrics as simply reflecting a growth of the effective number of degrees of freedom with 
temperature or energy scale - this has been suggested for the Dp-brane metrics in e.g. [36J. 
On the other hand, it has also been suggested that spin-glass phases of the random field Ising 
model can be governed by nontrivial hyperscaling violation exponent [3E]. Several groups 
have proposed different ways to model random or glassy phases with gravity or D-brane 
constructions [33 HHJ EHJ [50]; it would be very interesting if coarse-graining appropriately 



in any of these approaches, yielded metrics of the form (1.3) for reasons similar to those 
espoused in [3E]. 

Also, while here we followed an effective approach studying holography on slices at finite 
radius (associated to the cross-over scale />), one could imagine trying to take the limit rp — \ 
0. It would then be interesting to extend the methods of holographic renormalization [51] 
to metrics with hyperscaling violation. 

Finally, we uncovered here some novel "bottom up" holographic ground states with en- 
tanglement entropy intermediate between area law and extensive scaling, for d — 1 < 9 < d. 
It is well known that systems with extensive ground state entropy, like the duals of AdS% x M 2 
gravity theories (or, even simpler, theories of free decoupled spins), can yield extensive entan- 
glement entropy. It would be interesting to find candidate field-theoretic models which could 
yield the intermediate scalings we found here. In [52] and references therein, supersymmetric 
lattice models which are either "superfrustrated" (enjoying extensive ground-state entropy), 
or frustrated with large but sub-extensive ground state degeneracy, are described. It is quite 
possible that one can construct analogous lattice models with intermediate ground state 
entropies, giving rise to entanglement scaling like that of our new phases [53]. In fact, soon 
after this work was submitted, similar intermediate scalings of the entanglement entropy 
were found in field-theoretic models with impurities [53]. 
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A Metric properties 

This sections contains results on the Ricci and Einstein tensors for the general class of metrics 

dsl +2 = e 2 ^) (-e^f{r)de + ^ + dxf) . (A.l) 

The nonzero elements of the Ricci tensor are 
Rtt = \e 2B[r) f{r) ( [(d + 2)A\r) + 3B'(r)} f(r) + 2/(r) [(A'(r) + B> \r)){dA \r) + B'(r))+ 



+ 



A"(r) + B"(r)]+f"(r)) 



Rrr = -\f(r)- 1 ( [(d + 2)A'(r) + 3B'(r)} f (r) + 2/(r) [B'(r)(A'(r) + B'(r))+ 

+ (d+l)A"(r) + B"(r)] + f"(rf) 
% = (df(r)A'(r) 2 + f{r)A\r)B'{r) + A'(r)/'(r) + f(r)A"(rj) , (A.2) 
and the scalar curvature is 

i? = _ e -2A(r) ^ 2(rf + x)^^)/'^) + W'(r)f(r) + (d + l)/(r) [tL4'(r) 2 + 2A'(r)B'(r)] 

+ 2/(r) [^(r) 2 + (d + l)A"(r) + B"(r)] + /"(r)) . (A.3) 
As a result, the Einstein tensor Gmn = Rmn — \qmnR simplifies to 
G tt = -\de 2B ^f{r) (A(r)f'(r) + f(r) [(d - l)A(r) 2 + 2A"(r)} ) 
Grr = l -r\r)dA\r) (/(r) [(d + l)A'(r) + 2E'(r)] + /'(r)) 

^ = ^ " l)/(r)A'(r) 2 + 2d/(r)A'(r)B'(r) + 2f(r)B' (r) 2 + 2dA'(r)/'(r) + 

+ 33'(r)/V) + 2d/(r) A"(r) + 2f(r)B"(r) + /" (r)) . (A.4) 

We can impose the null energy condition, 

T^ U N»N V > 0, (A.5) 

where N^N^ = 0, on the Einstein equations G^ u = T^ u to derive constraints on the metric 
functions A(r),B(r), f(r). Choosing 

^ = ^77' N r = ^«»<P, N'=^ V , (A.6) 
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where 6 is an arbitrary constant, we get the constraints 

f(A' 2 + A'B' - A") > 

{dA' + 3B')f + 2f{dA'B' + B' 2 + 

from (p = and p> = | respectively. 



B") + /" > 



(A.7) 



B Massive propagators for general geodesies 

We consider the action 

S = -m [ dr r-^-^^^r-^-^T 2 + r 2 + x 2 



(b.i; 



where we have set A = r and r = it in equation (2.10), this time for both \Ax\ and |A 



nonzero. The integrated x, r equations of motion define two conserved momenta: 

r -(d-6)/d- 



y/l + x 2 + r- 2 ( z -Vf 2 '' 
r -{d-6)/d r -2{z-l)j. 



Vl + x 2 + r- 2 ( z - 1 )r 2 ' 
We can use these to rewrite x, f in terms of r and Hz, Il r : 



(B.2) 



and 



dx 
dr 

dr 
dr 



IV 



,-(*-!) 



' r -2{z-e/d) _ fJ2 r -2(z-l) _ TT2 



n T 



,(*-!) 



(B.3) 



(B.4) 



^ r -2(z-e/d) _ n2 r -2(z-i) _ n 2 .' 

Also, using the fact that at the turning point dr/dx, dr/dr = 0, we can derive a relationship 
between r t , IT^, I1 T , 

r -2( Z ^) _ n 2 r -2(,-l) _ n 2 = ^ (R5) 



Plugging equations (B.3) and (B.4) back into the action, we get an expression for the total 
geodesic distance, 

29/d-z-l 

^^^^=^^^^= (B.6) 

-2{z-6/d) _ Jj2 r -2(z-l) _ TJ 2 



S = —m / dr 



as a function of the conserved momenta. Equations (B.3), (B.4) and (B.6) can in principle 



be solved numerically for particular values of the critical exponents in order to find \Ax\, 
|At|, and S respectively in terms of 11^, I1 T and the turning point r t . These solutions, along 



with the constraint equation (B.5) can be used to find S(Ax, At). 
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